Influence of the mean flow on zonal flow generation 
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Excitation of zonal flow by the modulational instability in the presence of mean shear flow is 
considered. It is shown that the small amplitude mean flow favours the modulational instability, 
increasing the growth rate, whereas sufficiently strong mean shear significantly reduces the instability 
growth rate. 

PACS numbers: 



It is now widely recognized that zonal flows, i. e. az- 
imuthally symmetric band-like shear flows that depend 
only on the radial coordinate^ play a crucial role in 
regulating drift-wave turbulence and transport in toka- 
maks. It is now quite clear that zonal flows are gener- 
ated by modulational instability of drift waves^£&L£. 
It is important to distinguish the zonal flow from the 
mean shear flow associated with the mean radial elec- 
tric field. The latter can be driven and sustained in the 
absence of turbulence (by heating, fueling, momentum 
input, etc.), whereas zonal flows are driven exclusively 
by nonlinear wave interaction processes (via the mod- 
ulational instability) 1,2 . In contrast to smooth static 
mean flows, the zonal flow patterns have complex, pos- 
sibly random, spatial structure. It is well known that 
the presence of mean flow give rise not only to instability 
of the sheared layer (the Kelvin Helmholtz instability), 
but also to stabilization of other instabilities (ion temper- 
ature gradient driven modes, resistive interchange modes 
etc.) ^ 9 ' 10 ! 11 ' 12 . In this Brief Communication the influ- 
ence of the mean shear flow on excitation of zonal flow 
by the modulational instability is reported. 

We use a simple slab two-dimensional (x and y being 
the radial and poloidal coordinates respectively) model 
similar to the Hasegawa-Mima model, but taking into ac- 
count that the plasma density does not follow the Boltz- 
mann distribution for large scale motions^&i. We rep- 
resent the electrostatic potential as a sum of small scale 
fluctuating (p and mean (i^) quantities ip = ip-\-(ip), where 
(ip) = ip + cpo with ip being the zonal flow potential and 
Po corresponds to the mean shear flow. Equations de- 
scribing the coupled dynamics of large (zonal and mean 
flows) and small (drift waves) scale motions 
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(v) • VA^- v- VA(tp + ip ) = 0, (1) 



— A^+(v- VA^) =0, 



(2) 



where v = z x V<y5, (v) = vo(x)y + z x V<p and vq(x) 
is a profile of the mean flow, (...) denotes the averaging 
process. The normalizations of space in unit of p s , time 
in units of L n /c s , and potentials in units of T e p s /L n e 
are used in Eqs. (Q} and ([2]). Here notations are stan- 
dard, i.e. T e is the electron temperature, c s = ^jT e /m,i, 



p s = c s /VLi is the ion-sound Larmor radius, and L n is the 
characteristic length scale of the plasma inhomogeneity. 
In the linear approximation Eq. JT]) yields the ordinary 
dispersion relation for drift waves 
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The term in the angular brackets of Eq. ([2]) corresponds 
to the Reynolds stress induced by the small scale drift 
waves. 

The electrostatic potential <p of the drift waves is con- 
sidered as a superposition of a pump wave and two satel- 
lites, ip = ipd + <p+ + (p—. In what follows, we consider 
the case when the drift wave is propagating nearly in the 
poloidal direction, k x — and represent 

p> = p{x)e~ mt + v*(x)e mt , (4) 
tp d = $ oe l ( fe ^-^t) + $* e -i(fc„y-^kt) ) ( 5 ) 

<P± = ip±(x)e z{k ^- UJ±t) + ip* ± (x)e- i( - k v y - u±t \ (6) 

where u>± = u>k ± f2. Substituting Eqs. d4j) (J6j) into Eq. 
(|2|), we obtain 



nip = k y [ $0 
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(7) 



dx dx 

For the amplitudes of the satellites we have from Eq. ((T|) 
[fi + Lo k - k y v (x)} (l + kl - 
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where v'q = d 2 vo(x)/dx 2 . In the absence of mean shear 
flow (vq = 0), taking <p,ip+ ~ exp(iqx),we recover the 
previous dispersion relatio n 13 : 14 
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1 + k 2 + q 2 



(10) 
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FIG. 1: The growth rate versus k y for go = 0.05 and (a) mean 
flow amplitude V = and V = 0.04, drift wave amplitude 
| <3E> 1 = 0.02; (b) mean flow amplitude V = and V = 0.4, 
drift wave amplitude |$o| = 0.2. 

where 



(l + fe2)(i + fc 2 + g2) ' v ; 

which predicts the modulational instability (f2 2 < 0) if 
the drift wave amplitude |$o| is above the correspond- 
ing threshold. Under q <C k y maximum growth rate 
is achieved at q opt = k y \$ \(l + fc 2 ) 3/2 and j max = 
^|<I> | 2 (l + fc 2 ). 

In the presence of mean flow (vq ^ 0) equations (JZ))-© 
can be rewritten as the generalized eigenvalue problem 

A -A\ (l 0\ 

B C $ = fi 0£ * (12) 
—B D J \Q E) 

where $ = [ip, ip + ,(p*_), and 

A = -M>o^ B = -^o{k^ + ^), (13) 

C = ky + k y v'o - (uj k - k y v )E, (14) 
D = (wk + k y v )E - k y - kyVQ , (15) 

We assume that the profile of the mean shear flow has the 
form vq(x) — V tanh(qo^). Employing a finite differenc- 
ing approximation, we numerically solved the eigenvalue 
problem (fl"2"|) . The instability corresponds to eigenvalues 
with ImQ ^ 0. Figures [TJa) and (b) show the maximum 
growth rate 7 = max|Im£7| of the modulational insta- 
bility modified by the mean shear flow as a function of 



the poloidal wave number k y . Dashed lines correspond to 
the case of the absence of mean flow. It is seen that the 
small amplitude mean flow [Fig. [lja)] favours the modu- 
lational instability and zonal flow generation, increasing 
the growth rate, whereas sufficiently strong mean shear 
[Fig. [3Jb)] significantly reduces the instability growth 
rate. In Fig. [5] we plot the growth rate for fixed k y and 
several values of the drift wave amplitude $0 1 a s a func- 
tion of the mean flow amplitude V. The growth rate 7 
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FIG. 2: The growth rate versus the mean flow amplitude V 
for three different values of the drift wave amplitude: |<&o| = 
0.02, 0.05, 0.1. Here k y = 1 and go = 0.05. 



initially increases with increasing V, and thus the pres- 
ence of mean flow has a destabilizing effect. Then, 7 
attains its maximum at some V (V ~ 0.08 for k y = 1 
and all |$o|)- This corresponds to the most unstable 
case. Then, the growth rate decreases as V increases, 
and when the mean flow amplitude exceeds some critical 
value, the growth rate becomes smaller then in the case 
of the absence of mean flow. Further increasing of V 
can significantly reduce the growth rate. Thus the pres- 
ence of mean flow with sufficiently strong amplitude has 
a stabilizing effect on generation of zonal flows by the 
modulational instability. 
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